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ABSTRACT 

We present results from Chandra observations of the z\ = 0.58 cluster lens 
SDSS J1029+2623, which is presently the largest separation gravitationally 
lensed quasar. We clearly detect X-ray emission both from the lensing clus- 
ter and the three lensed quasar images. The cluster has an X-ray temperature of 
kT = 8.ll^2 keV and luminosity Lx = 9.6 x 10 44 ergs -1 . Its surface brightness 
is centered near galaxy G2 and it is elongated East- West. We identify a subpeak 
North- West of the main peak, which is suggestive of an ongoing merger. Even 
so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears 
to be consistent with the Einstein radius of the system. We find significant ab- 
sorption in the soft X-ray spectrum of quasar image C, which can be caused by 
an intervening material at either the lens or source redshift. The X-ray flux ra- 
tios between the quasar images after correcting for absorption are in reasonable 
agreement with those at optical and radio wavelengths. This implies that mi- 
crolensing effect is not significant for this system and dark matter substructure 
is mainly responsible for the anomalous flux ratios. 

Subject headings: galaxies: clusters: general — gravitational lensing: strong - 
quasars: individual (SDSS J1029+2623) — X-rays: galaxies 
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Introduction 



The gravitational lens SDSS J1029+2623 consists of three images o 



quasa r produced by a foreground cluster at zi ~ 0.58 fllnada et al 



2006 



8b Zc 



2.197 



Oguri et al 



20081 ). With a maximum image separation of 22'.'5, this is the largest lens among ~ 130 
gravitationally lensed quasars known to date. The lensing interpretation is secure from 
multiple observational facts, including the remarkably similar optical spectra with mini 
broad absorption lin e (BAL) features in their emission lines, similar radio loudnesses 



( IKratzer et al. 



with weak lensing ( jOguri et al 



201 if), and the identif ication of the lensing cluster in the optical and 



20121 ). Furthermore, this lens system represents a rare 



example of a naked cusp lens. In this configuration, only three images with similar 
brightnesses are produced on the same side of th e lens potential, which has been predi cted 



to be fairly common amon g cluster-scale lenses ( jOguri fc Keeton 



2004 



Li et al. 



2007 



Minor fc Kaplinghat 



2008|). 



SDSS J1029+2623 is one of two examples of large-separation quasar lenses produced 



by massive c 



( jlnada et al. 



usters of galaxies. The other large-separa t ion lens is SPSS J1004+4112 



2003 



Oguri et al. 



2004 



Sharon et al. 



2005 



Inadaetal 



20051 ). In fact, these 



clusters can be regarded as unique examples of strong-lens-selected clusters in the sense 
that they have been identified by searching for strong lenses in a large sample of quasars. 
Therefore, detailed studies of the lensing clusters may provide useful insights into the nature 
and diversity of intermediate redshift clusters. In particular, (1) do the strong- lens-selected 
clusters follow empirical scaling relations derived from optical/X-ray selected clusters? and 
(2) what is the dynamical state of these clusters? 

Chandra and XMM-Newton observations of the other large-separation quasar lens, 
SDSS J1004+4112, have shown that the lensing cluster is relaxed, with the X-ray 
temperature and luminosity being consistent with the predictions of empirical scaling 
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relations within the observed scatter (jOta et al.l l2006bt lLamer et al.l 120061 ). This 
interpretation is also suppor ted by the r ecent comparison of the X-ray measurements with 



detailed strong lens models ( IQguri 



2010|). There is excellent agreement on the centroids and 
position angles between the inferred dark matter distribution, the X-ray surface brightness, 
and the brightest cluster galaxy. 

The available observational data suggest that the lensing cluster of SDSS J1029+2623 
may be quite complicated. For instance, the existence of two bright elliptical galaxies at the 
center with a large velocity difference (~ 2800 kms -1 ) between them indicates a possible 
merger. High-resolution Chandra imaging and spectroscopy are essential for studying the 
dynamical state of the lensing cluster, based on detailed comparisons between the gas 
distribution probed by the X- ray measurements and the mass distribution obtained from 



strong and weak lensing (e.g. 



Clowe et al 



2006|) 



In addition to the properties of the lensing cluster, the high-resolution Chandra 
observations provide an additional probe of the anomalous flux ratio between the quasar 
images. In particular, the large flu x difference between the close quasar image pair (images 



B and C) is observed in the radio (IKratzer et al.l 1201 if ), which implies that a large part of 



ensing clust er. 



the flux anomaly is caused by (sub-)structure in the mass distribution of the 
Howe ver, there is a clear wavelength dependence to the optical flux ratios ( IQguri et al.l 
20081 ). which indicates that chromatic microlensing or dust extinction must also affect the 
third quasar image. 

In this paper, we present results from Chandra observations of this unique cluster-scale 
quasar lens. We successfully detected X-ray emission from both the lensing cluster and the 
three quasar images, which are discussed in §[3] and respectively. We compare our results 
with observations at other wavelengths, as well as strong lensing information from the 
quasar images. Throughout the paper we adopt a cosmological model with matter density 
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Qm — 0.27, cosmological constant Cl\ = 0.73, and Hubble constant Hq = 70 kms _1 Mpc _1 . 
At the cluster redshift of z — 0.584, 1" corresponds to 6.69 kpc. Unless otherwise specified, 
quoted errors indicate the 90% confidence range. 



2. Observations 



SDSS J1029+26 23 was observed for 60 ks with the Chandra Advanced CCD Imaging 



Spectrometer (ACIS; 



Gar mire et al 



[20031) on 2010 March 11 (Observation ID: 11755). The 
data were obtained with the ACIS S3 CCD operating in VFAINT mode. This CCD has a 
1024 x 1024 pixel format with an image scale of / /492 pixel -1 . The target was offset from 
the nominal aim point with a Y-offset of —1', which has little effect on the spatial resolution. 
The CCD temperature during the observations was — 120°C. The data were processed 
using the standard software packages CIA0 4.3 and CALDB 4.4.1. The background was 
stable during the observation and the net exposure time after applying standard light-curve 
screening cuts was 55815 s. 

The ACIS image in the 0.5-7 keV band is shown in Figure [TJ The extended cluster 
emission and the three images of SDSS J1029+2623 are all detected. The astrometry offset 
is found to be negligibly small (0'.'2) from the comparison of quasar A position with the 
optical image. The cluster emission exhibits a somewhat elongated feature, which we will 
study in detail in §3.21 In the ACIS-S3 field, 42 point-like sources, including the three 
quasar images, were detected with the wavdetect algorithm for a significance threshold 
parameter of 10~ 6 . The X-ray centroid position (10:29:12.47, +26:23:33.2), measured in 
X2l agrees well with the position of the central galaxy G2 at (10:29:12.48, +26:23:32.0; 



Inada et al. 



20061 ). 
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3. Lensing cluster 



3.1. Spectral analysis 



We derive the spectrum of the cluster component by extracting the data from a 
circular region with a radius of 1' centered on the X-ray centroid. We exclude 5" radius 
regions around each quasar image. The background was estimated in a surrounding annulus 
(2f2 < r < 2f5). The source count rate within r < V is 0.068 ± 0.001 counts s" 1 . We fit 
the cluster spectru m in the 0.5-7 keV band with the APEC thin-thermal plasma model 



( jSmith et al. 



20011 ) utilizing XSPEC version 12.6.0q. The Galac tic hydrogen co 



was fixed to Nr = 1.67 x 10 20 cm 2 based on the LAB survey (IKalberla et al. 



umn density 



2005|). Given 



the good agreement between the X-ray emission centroid an d the position of G 2 galaxy (see 



below), we set the cluster redshift to that of G2 (z = 0.584; 



Oguri et al. 



20081 ) 



The spectrum is shown in Figure |2j From the spectral fits the X-ray temperature is 
constrained to be kT = 8.ll 2 '2 keV and the metal abundance to be Z = 0.21Z Q (the 90% 
upper limit is O.44Z ). This high temperature indicates that the lensing cluster is indeed 
massive enough to be capable of producing the large-separation lensed quasar images. The 
reduced chi-square of the best-fit model is x 2 /dof = 138/139, and the parameters of the fit 
are summarized in Table [TJ 

The absorption-corrected 0.5-7 keV flux is 4.6 x 10 _13 ergs _1 (r < 1'). The bolometric 
X-ray luminosity within r 50 o is estimated as Lx = 9.6 x 10 44 ergs -1 , where r 500 is defined 
as the radius within which the average matter density is equal to A c = 500 times the 
critical density of the Universe at the cluster redshift. The 0- model analysis in § 13.21 
yields r 50 o = 1.09 Mpc. This luminosity is lower than the mean value expect ed from th e 



lumino sity-temp er at ure relation of distant clusters, Lx = 4.0 x 10 45 ergs 1 (lOta et al. 



2006a)), but it is within the observed scatter of the Lx — T relation. 
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To investigate the radial temperature profile, we also analyze spectra for the inner 
r < 0f25(~ 100 kpc) region and the outer 0'25 < r < V region of the lensing cluster. 
The fitted temperatures are kT = Q.Stl'l keV and 6.8l];£ keV, respectively. While the 
temperature profile is consistent with being constant with radius, the possible decrease of 
the temperature toward large radii implies that t he lensing cluster la cks the radiatively 



cooling core commonly seen in unrelaxed clusters ( lArnaud et al. 



2010h . 



3.2. Image analysis 



We first fit a single one-dimensional radial surface-brightness profiles to the extended 
X-ray emission. We centered the models at the X-ray centroid and azimuthally averaged 
the 0.5-5 keV image. We corrected for vignetting and the detector response using the 
exposure map calculated for the spectral energy distribution of the cluster. The image was 
rebinned by a factor of two so that the pixel scale is 0'/98 (6.58 kpc at the cluster redshift). 
The quasar images were masked out to 5" when calculating the profile. 

We fit the radial profile with the following two models: (1) the conventional /3-model 
S(r) = Sq[1 + (r /r r ) 2 ]~ 3 ^ +1 / 2 , and (2) the profile derived from the u niversal mass pro file 



proposed by 



Navarro et al. 



(119971 ) plus isothermality of the cluster (ISuto et al 



1998 



hereafter NFW-SSM). The background level, which was assumed to be a constant, is 
simultaneously fit. As illustrated in Figure [31 we find that both models can fit the observed 
radial profile reasonably well, with reduced chi-square values of % 2 /dof = 183.2/196 and 
178.9/196 for the /3-model and NFW-SSM model, respectively. The fitted values of the 
/3-model parameters are f3 = 0. 721^04 an d r c = 16.41];;® arcsec(= 109l] 2 kpc), while the 
NFW-SSM model parameters are B = 12.21] J and r s = 78l]^ arcsec(= 519l^ 4 kpc). 

In the /3-model, the central electron density is n e o = (1.5 ±0.1) x 10" 2 cm~ 3 . This 
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implies a radiative cooling time at the cluster center of t coo i = 5.8^q 9 Gyr, which is shorter 
than the Hubble time but comparable to the age of the Universe at the cluster redshift, 
^age = 8.2 Gyr. Radiative cooling should not be important in the cluster core region. 

The Chandra image suggests that the cluster emission is somewhat elongated along 
the East-West direction. Moreover, there is another emission peak about 10" North- West 
of the cluster centroid (Fig. [I]) . These motivate us to explore the X-ray surface brightness 
with two-dimensional models, adopting both the single- and double-component elliptical 

i ~\ is defined by 



/9-models. The elliptical /3-mode 



S(x,y) = S 



1+1 - 



2 

r 



-3/3+1/2 



C, (1) 



r(x,y) = [{(x - x )cos9 + (y - y ) sin9} 2 (l - e) 2 

+{(y ~ Vo) cos 6-{x- x Q ) sin 6} 2 ] 1/2 /(l - e), 

where r c is the core radius, x and y define the cluster center, e is the ellipticity, 
9 is the position angle, So is the amplitude at the center. The background level 
per image pixel was fixed to the value obtained from the one-dimensional analysis, 
C = 3.3 x 1CT 9 counts s -1 cm -2 . 

We first fit a single-component elliptical /3-model. The maximum-likelihood fit was 
performed with Sherpa, where the exposure map was included to convolve the model image 
with the telescope and detector responses. We find that the fit residuals show significant 
excess emission in the North- West region, as shown in Fig. HI The parameters and results 
for the single- component model are summarized in Table [3j 

We next consider a two-component model in order to determine the X-ray emission 
profile of the ICM more precisely. The model consists of two elliptical /3-models, one 



x beta2d model in the Sherpa fitting package 
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centered on the cluster and the other on the North- West component, plus the constant 
background. For the North- West component, the parameter /3 2 and position angle 82 were 
fixed to 0.7 and 0, respectively, due to the limited photon statistics. The results from 
fitting the central 3'3 x 3'.3 region are shown in Table [3] and Figure [51 In order to check the 
goodness of the fit, we re-binned the image into two single-dimensional profiles in the Right 
Ascension and Declination directions (Fig. [5J bottom right) and calculated the x 2 values 
between the model and data profiles within the central ±1' region to find that they are 
sufficiently small, x 2 /dof = 131/122 and 148/122 for the x- and y-directions, respectively. 
The best-fit cluster-center position of (10:2 9:12.47, +26:23:33.2) is consistent with G2 



10:29:12.48, +26:23:32.0; 



Inada et al. 



20061 ) within the errors. The ellipticity of the main 
cluster component is measured to be e = 0.214 ^°'°^ (axis ratio q = 0.786), which is quite 



Hashimoto et al. 



2007 



Kawaharall2010l : 



ticities ( 


De Fili 


DDis et al. 


2005 


Lau et al. 


2012 


)• 



Flores et al. 



2007 



Note that there appears to be no optical counterpart for the NW peak at (10:29:11.9, 
+26:23:37.4). For the APEC model, the X-ray spectrum extracted from r < 5" around 
the NW peak with subtracted background from the surrounding 5" < r < 10" region 
yields kT = 3.8(> 1.7) keV and an unabsorbed bolometric luminosity of 1.5 x 10 43 ergs -1 . 
Adopting the estimated bolometric luminosity, a temperature of 1.1 ke V is obtained from 



2006ah . 



the luminosity-temperature relation for such small core radius clusters (lOta et al. 
which is significantly lower than the temperature measured from the Chandra observation. 
This might be the signature that the NW component is in the process of merging with the 
main cluster or has recently undergone a merger. These parameters are reminiscent of t he 



20041 ). 



emission peak detected in the line-of-sight merging system CL0024+17 (lOta et al. 
In addition, there is a possible jump in the X-ray surface brightness distribution near the 
NW peak (Fig. ED) . The s e observed feature s suggest that the NW peak may be a cold front 



(IMarkevitch et al. 



2000 



Owersetal 



20091 ). although the limited photon statistics prevent 
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a more detailed exploration. 



3.3. Cluster mass distribution 



Under the assumption of hydrostatic equilibrium and spherical symmetry, we can 
infer the mass distribution of the lensing cluster from the X-ray temperature and 
surface-brightness profiles. We assume isothermality because we did not detect a significant 
radial dependence of the X-ray temperature. Figure M shows the cylindrical cluster mass 
projected within a radius r derived from each surface mass distribution profile. We find 
that both models yield consistent mass profiles within the measurement errors. 



For the NFW-SSM model and the critical overdensity A c = A^ 



187T 2 fi°- 427 



(INakamura fc Suto 



19971 ). the virial mass and the concentration parameter are constrained 



to be M v ir 



1.57+°™ x 1O 15 M and c vir = 4.52if!«. However, these values should be 



-0.39 ~ ±u J "0 allLt u vir — ^-^-0.58- 

taken with caution because they rely on the extrapolation of the profiles well beyond the 
radii where X-ray properties are measured. For comparison with other mass measurements, 
it may be better to use masses defined for larger critical overdensities A c . In Table H] we 
summarize the mass and radius estimated for the virial overdensity and overdensities of 
A c = 500, 1000, and 2500, for both the /3-model and the NFW-SSM model. 

Next we compare the mass distribution derived from the X-ray analysis with those 
from strong lens models. An issue for strong lens modeling for this cluster has been the 
uncertainty of the mass centroid, be cause the naked cusp configuration poorly constrains 



the center of the mass distribution ( Inada et al. 



2006 



Oguri et al. 



20081 ). Given that the 



measurement of the X-ray centroid agrees well with the location of one of the brightest 
cluster galaxies (galaxy G2), it now seems safe to assume that the centroid coincides with 
galaxy G2. We fit the positions of the three quasar images assuming an NFW profile plus 
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201(1 ) . We find that the best-fit 



external shear model using the software glafic (lOguril 
model predicts the Einstein radius of 8^ = 18'/3 for the source redshift of z s = 2.197. This 
translates into a cylindrical cluster mass of M 2 d(< 122kpc) = 9.30 x 1O 13 M , which is in 



good agreement with the enclosed mass derived from the X-ray data (see Figure [6]). 



4. Quasar images 

4.1. X-ray properties 

Given the small 1'.'8 image separation of images B and C (see Figure [7]), we measure 
the X-ray fluxes and spectra of the quasar images using a circular aperture with a radius 
of 0'.'9 for all the quasar images. The observed counts in the 0.5-7 keV band are 687, 1073, 
and 295, for images A-C, respectively. The net contamination from the background and the 
cluster is estimated to be only a few percent. Similarly the contamination from image B in 
the extraction radius for image C is estimated to be ~ 10 counts, which is smaller than the 
la statistical error of the flux of image C. We do not find any clear sign of time variability 
for light curves produced with time resolutions of 2048 or 4096 s. 

The combined spectrum of the three lensed images A, B, and C is shown in the left 
panel of Figure El We fit the spectrum with a power-law model. The Galactic absorption 
is again fixed at Nr = 1.67 x 10 20 cm -2 . We do not find significant intrinsic absorption 
in the total spectrum. The power-law model is acceptable at the 90% confidence level 
(% 2 /dof = 57.1/63) and the power-law index of F = 1.55 ± 0.06 agrees well with the mean 



value for mini-BAL quasars observed with Chandra (jGibson et al.ll2009l ). The corresponding 
total luminosity is 5.3 x 10 45 ergs s _1 in the 2-10 keV band. We do not detect any significant 
emission from the neutral iron Ka line at 2.0 keV in the observed frame (or 6.4 keV in the 
rest frame). The 90% upper limit on the equivalent width is EW < 170 eV in the quasar 
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rest frame for a power-law combined with a Gaussian line model at 2 keV. 

Next we study the spectra of the individual quasar images. We first check the hardness 
ratios defined by the count rate ratio in the 2-7 keV and 0.5-2 keV bands. The result 
indicates that a spectral shape of image C differs from A and B because the hardness ratio 
for image C (0.48 ± 0.06) differs from those of images A (0.37 ± 0.03) and B (0.36 ± 0.02). 
Indeed, the spectra of images A-C plotted in the right panel of Figure [8] clearly indicate 
that the emission of image C is heavily absorbed in the soft band. To quantify this, we fit 
the spectrum with a power-law model absorbed by intervening cold material in addition to 
the Galactic interstellar medium. Since the location of the cold absorber is not known, we 
consider cases with the absorber located at either the quasar redshift of z s = 2.197 or at the 
lens redshift of z\ = 0.584. 

The result of the power-law plus the absorption model fits are summarized in Table [51 
While the inclusion of the additional absorption does not improve the best-fit x 2 value for 
images A and B, the fits to image C are significantly improved by including the absorption. 
We find that the improvement is significant at the > 99.98% confidence level according to 
the F-test (Ax 2 = 11 for one additional parameter), although the redshift of the absorber 
is not constrained. We also tried a simultaneous fit to the three spectra with a common 
intrinsic spectral index T, because the intrinsic photon indices are statistically consistent 
among the three images. The result shown in Table [6] again indicates that the extra 
absorption is significant only for image C. 

4.2. Comparison with optical and radio flux ratios 



In Table U\ we compare flux ratios between quasar imag es measured in the C handra 
X-ray images with radio and optical flux ratios presented in iKratzer et al.l ( 1201 ll ). We 
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find that the absorption corrected flux ratio C/B in X-ray is larger than the optical flux 
ratios and is close to the radio flux ratio. We note that the optical C/B flux ratio becomes 
C/B ~ 0.4 once the dust extinction is corrected (see below), which is also consistent with 
the radio and absorption-corrected X-ray flux ratios. The nearly consistent flux ratios 
between the radio, optical, and X-ray observations imply that microlensing effect is not 
significant for this system, particularly because X-ray flux ratios are mo re easily affected by 



microlensing (e.g., 



Ota et al. 



2006b 



Pooley et al. 



2007 



Dai et al 



2010|). 



The origin of the absorption may be inferred from the dust-to-gas ratio. Based on 
the optical spectra of images B and C in 



Oguri et al 



(120081 ). we estimated a color excess 



of E(B - V) ~ 0.07 for absorption at z s = 2.197 and E(B - V) ~ 0.17 for absorption 
at zi = 0.584. These imply dust-to-gas ratios of E(B — V)/Nh ~ 2 x 10 _24 magcm 2 and 
E(B — V)/Nh — 3 x 10 _23 magcm 2 at the source and lens redshifts, respectively. The 



higher value found for dust in 



lens systems (IDai &: Kochanek 



;he le ns is broadly consistent with measurements for other 



20091 ). while the value for dust at the source redshift is much 



lower than expected from this observed correlation. However, suc h low dust 



;o-ga s ratios 



20061 ). which 



appear to be common for intrinsic absorption features of quasars (lHall et al. 
can be explained by dust evaporation due to irradiation by the central source. Thus we 
cannot draw any firm conclusion about the origin of the absorption fro m this analysis. Both 



the mini-BAL features in the quasar spectrum (see 



Oguri et al 



of dust in the intracluster medium particularly near the center (IChelouche et al 



Muller et al. 



2008 



Bovy et al. 



2008 



2008 ) and the deficienc y 



2007 



Kitayama et al. 



20091 ) suggest that intrinsic absorption 



at the quasar redshift is the more likely scenario, although the anomalous flux ratios imply 
a faint galaxy near image C which might host the dust at the lens redshift. If the former 
interpretation is correct, the different absorptions between image B and C probe spatial 



variations o 



(I Green 



20061 ) 



the absorption content, possibly related to the structure of disk wind outflows 
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5. Summary and Discussions 

We have presented results from Chandra observations of SDSS J1029+2623, the 
largest-separation gravitational lens system, produced by a cluster of galaxies at z — 0.58. 
We detect significant X-ray emission from the lensing cluster as well as three lensed quasar 
images. 

From the extended cluster emission, the X-ray temperature and bolometric luminosity 
are constrained to be kT = 8.lli2 keV and Lx = 9.6 x 10 44 ergs -1 , which is consistent with 
the luminosity-temperature relation of distant clusters. While the estimated cooling time 
is shorter than the Hubble time, the lack of a significant temperature drop at the cluster 
center suggests that radiative cooling is not important in this system. From the image 
analysis, we have found that the cluster emission is elongated along the East- West direction 
and has a subpeak in the North- West region. This indicates that the system has undergone 
a merger. However, the modest ellipticity of the gas distribution found in the main cluster 
component indicates that the irregularity of mass distribution may not be large. We have 
reconstructed the mass profile of the lensing cluster, assuming isothermality and hydrostatic 
equilibrium to find that the X-ray mass within the Einstein radius or ~ 120 kpc is in a good 
agreement with that expected from the strong lensing. The result further supports the idea 
that there is not any significant departure from the equilibrium state in the cluster core. 

The X-ray spectra of the three quasar images above 2 keV are well represented by a 
power-law with a photon index of F ~ 1.55 that is typical of mini-BAL quasars. However, 
significant soft X-ray absorption due to the cold material is detected only in image C. After 
correcting for the X-ray absorption, the X-ray flux ratios between the quasar images are 
roughly consistent with those in optical and radio, which suggests that microlensing is not 
significant for this system. This indicates that the flux ratio anomaly s een in this lens 



system has to be due to substructure (see, e.g., 



Kochanek fc Dalai 



20041 ). It is difficult to 
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distinguish the possibilities between absorbers at the quasar redshift and at the cluster 
redshifts from the X-ray spectrum or the dust-to-gas ratio. 

Both the surface brightness profile of the lensing cluster and the anomalous flux 
ratio are better understood by comparing the result with accurate mass modeling results 
obtained from gravitational lensing analysis. For this purpose we have recently obtained 
Hubble Space Telescope images of this system. Lensing analysis results of the HST data and 
a comparison with the Chandra results will be presented in a forthcoming paper. 
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Table 1: APEC Models of the Overall Spectrum 



Parameter Value (90% error) 

N u (cm- 2 ) 1.67 x 10 20 (F) 

kT (keV) 8.1 (6.8 - 10.1) 

Abundance (Z Q ) 0.22(< 0.44) 

Redshift 0.584 (F) 

Normalization 51 7.31 (6.93 - 7.68) x 10- 4 

X 2 /dof 47.5/45 
/x.o.5-7 (ergs-^m- 2 ) 13 4.6 x lO" 13 

Lx,o. 5 -7 (ergs- 1 ) c 5.5 x 10 44 

^x,boi (ergs" 1 ) d 9.6 xlO 44 



a Normalization factor for the APEC model, / n e nndV/4:irD A (l + z) 2 in 10 14 cm 5 , where Da is angular size 
distance to the cluster. 

b Galactic absorption-corrected X-ray flux in the 0.5-7 keV band within r < V . 
c Galactic absorption-corrected X-ray luminosity in the 0.5-7 keV band within r < 1'. 
d Bolometric luminosity within r 500 . 
(F) Fixed parameters. 
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Table 2. Results for Circularly Symmetric Models 



Model 



Parameter 



Value 



/9-model 



S (counts s- 1 cm" 2 kpc- 2 ) a 3.57^ x 10" 



r c (arcsec/kpc) 
X 2 /dof 



n 7 9 +0.05 
u - ' z -0.04 

16.4+H/1091} 2 
183.2/196 



NFW-SSM model S (counts s" 1 cm' 2 kpc~ 2 ) a 4.03lg;|J x 1(T 9 

B 12.2±\% 

r s (arcsec/kpc) 78±^/519i^! 4 
X 2 /dof 178.9/196 
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Table 3. Results for Elliptical Models 



Model Parameter Value 

Single elliptical /3-model Zo,i,S/o,i 10:29:12.395, +26:23:33. 71 a 

S (counts s- 1 cm" 2 kpc" 2 ) 3.57+8:1 x 1(T 9 

/3 0.70±g:8S 



r c (arcsec/kpc) 17.6+^/118+^ 

0.214^ 

e 0.3618;" 



Double elliptical /3-model »o,i,yo,i 10:29:12.471, +26:23:33.24 b 

S ,i (counts s+ cm" 2 kpc" 2 ) 3. 15±g;|? x 10~ 9 

0.71+ - 04 



Pi u./±_ 004 
r c ,i (arcsec/kpc) 19.0 +2 ;<j/127 + ^ 

0.21418;8^ 



0! 0.36+ - 09 



-0.035 

-0.09 
'-0.09 

x ,2,yo,2 10:29:11.853, +26:23:37.37 c 

S , 2 (counts s- 1 cm" 2 kpc" 2 ) 3. lt 2 ;^ x 10~ 9 

0.70(F) 

r Cj2 (arcsec/kpc) 3.9+^/26 +2 j 
e 2 0.19(< 0.35) 

#2 0.0(F) 
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a The 90% errors are il'.'l for x and ±1"0 for y . 
b The 90% errors are ±l'/4 for x ,i and il'.'l for y 0)1 . 
c The 90% errors are ±l'/7 for x ,2 and ±1'.'5 for y ,2- 
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Table 4: Masses for Different Critical Overdensities 



Model 


M ■ a 

1 VI vlr 


t- . b 
' vir 


M 500 a 


r 500 b 


^iooo a 


r iooo b 


^2500 a 


r 2500 b 


NFW-SSM 


12.951*;°? 

15.07^3.93 


2 i n +0 - 24 

9 O/ i+0.28 
z - z ^-0.21 


6.681?;" 
9.28+^f 


i.uy_ 09 

1 99+O.I5 
i - zz -0.12 


4-65±i;f? 
6.45l?$ 


n 77 +0.09 
u - ' ' -0.07 

U -^ D -0.08 


9 a-t +1.13 
z " o± -0.68 

o r +1.62 
°- JZ -0.99 


o.48i8:S 

0.521™ 7 



Q In units of 10 14 M ( 
b In units of Mpc. 
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Table 5: Absorbed Power-law Models for the Quasar Images 





No extra absorption 


Absorption at z — 


2.197 


Absorption at z = 


0.584 


Name 


r 


X 2 /dof 


r 




X 2 /dof 


r 


iV H a 


x 2 /dof 


A 


1-5618:1? 


31.3/33 


i-56±8:i? 


0.00(< 0.70) 


31.3/32 


i.56t8:i? 


0.00(< 0.12) 


31.3/32 


B 


i-68±8:8g 


23.2/33 


1-6818:58 


0.00(< 0.33) 


23.2/32 


1-6818:58 


0.00(< 0.06) 


23.2/32 


C 


i O1+0.14 
i - zi -0.14 


25.1/25 


1 68 +a32 

l.UO_ 28 


q I q+2.24 
°- 10 -1.68 


13.5/24 


1 71+0.34 

L - ' x -0.30 


54+ - 40 

u -°^-0.30 


14.1/24 


Total 


i prcr+0.06 
i - od -0.06 


57.1/63 


i cc+0.10 
i - do -0.06 


0.00(< 0.46) 


57.1/62 


1 cr+O.lO 
i -°°-0.06 


0.00(< 0.08) 


57.1/62 



a The hydrogen column density in the unit of [10 22 cm 2 j. 



Table 6: 


Joint fits to the Quasar 


Spectra 








Absorption at z = 


2.197 


Absorption at z — 


0.584 


Name 


r iv H a 


X 2 /dof 


r iv H a 


xVdof 


A 
B 

C 


1.6518:8s 0.24(< 0.89) 

0.00(< 0.25) 

2 Q7+ L38 
z - y ' -1.08 


69.4/90 


1.6518:8s 0.04(< 0.15) 
0.00(< 0.04) 

o.49l8:^ 


70.0/90 



a The hydrogen column density in the unit of [10 22 cm 2 ]. 
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Table 7. X-ray and optical properties of SDSS J1029+2623 



Name 



Fx a,( 



hardness 



/(X) b < d f(opt,V) d /(opt,/) d /(radio) d 



A 0.90 ±0.03 (0.89) 0.37 ±0.03 -1.32 0.66 (0.64) 0.96 0.95 0.80 ± 0.06 

B 1.37 ±0.04 (1.34) 0.36 ± 0.02 -1.26 = 1 = 1 = 1 =1 

C 0.49 ±0.03 (0.40) 0.48 ± 0.06 -1.14 0.36 (0.27) 0.16 e 0.24 c 0.46 ± 0.05 



a Absorption-corrected, 0.5-7 keV X-ray flux in units of 10 _13 erg s _1 cm -2 and the l-er error. 
b The value in the parenthesis are those without absorption corrections. 

c Hardness ratios defined by the count rate ratios in 2 — 7 keV and 0.5 — 2 keV bands. Errors indicate 
68% confidence limits. 



J Flux normalized by the flux of image B. The optical and ra dio flux ratios are taken from 



Kratzer et al 



(j2011r ). Statistical errors on the optical flux ratios are < 2% (jOguri et al 



2008) 



°The optical C/B flux ratio becomes C/B ~ 0.4 once the dust extinction is corrected. 



-26 - 




Fig. 1.— Adaptively-smoothed ACIS-S3 image of SDSS J1029+2623 in the 0.5-7 keV band. 
The A-C quasar images of SDSS J1029+2623 and the extended emission from the lensing 
cluster are seen. The X-ray centroid of the cluster component and the North- West subpeak 
are marked with an "x" while the positions of the two most luminous galaxies, Gl and G2, 
are marked with a "+". The close B/C image pair is resolved in the raw, unsmoothed image 
(see §0| . 
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Fig. 2.— Chandra ACIS-S3 spectrum of SDSS J1029+2623 (r < 1') fitted with the APEC 
model. In the top panel, the crosses denote the observed spectrum and the histogram function 
shows the best-fit model convolved with the telescope and detector response functions. The 
bottom panel shows the residuals of the fit in units of a. 
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Fig. 3. — Results of the X-ray surface brightness profile fits with (a) the /3-model and (b) 
the NFW-SSM model. In each panel, the crosses show the observed surface brightness in 
the 0.5-5 keV band and the solid line shows the best-fit model. The background is shown 
with the horizontal dotted line. The vertical dashed line shows the observed extent of the 
diffuse X-ray emission, r x (see text). The bottom panel shows the residuals in units of the 
local noise. 
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Fig. 4. — Residuals from fits with a single elliptical /3-model. The image is smoothed with a 
a = 2" Gaussian, and excess emission over the best-fit model is shown by the black contours 
with linear spacing. 
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Fig. 5. — Results for the fits with the double elliptical /3-model. Top left: Chandra ACIS-S3 
image of the central 3'.3 x 3'3 region of SDSS J1029+2623 in the 0.5-5 keV energy range. 
Top right: The best-fitting elliptical /3-model, overlaid with logarithmically spaced intensity 
contours. Bottom left: Residuals from the fit. Bottom right: Residuals along East- West 
(top) and North-South (bottom) lines through the cluster center in units of significance. 



31 



t 1 1 — i — i — i i i | 1 1 1 — i — i — i i i | r 




0.01 0.1 1 

r [h^Mpc] 

Fig. 6. — Enclosed mass of the lensing cluster, Mx, for the /3-model (black) and the NFW- 
SSM model (red). The dotted lines indicate the 90% error ranges. Note that Mx is a 
cylindrical cluster mass projected within a radius r. The mass derived from gravitational 
lensing is also shown for comparison. The meaning of the vertical dashed line is the same as 
Figure El 




Fig. 7.— Raw ACIS-S3 image of SDSS J1029+2623 in the 0.5-7 keV band. The spectral 
regions for quasar images A-C are indicated with the circles. The meaning of the labels are 
the same as Figure [H 
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Fig. 8. — Left: total quasar spectrum fit with the power-law model (solid line). Right: 
Spectra of lensed quasar images A (black), B (red), and C (green) fit with a common intrinsic 
slope T. 



